Abstract Objective: Cleft palate (CP) is a congenital birth defect caused by the failure of palatal fusion. Little is known about the potential role of DNA methylation in the pathogenesis of CP. This study aimed to explore the potential role of DNA methylation in the mechanism of CP. Methodology: We established an all-trans retinoic acid (ATRA)-induced CP model in C57BL/6J mice and used methylation-dependent restriction enzymes (MethylRAD, FspEI) combined with high-throughput sequencing (HiSeq X Ten) to compare genome-wide DNA methylation profiles of embryonic mouse palatal tissues, between embryos from ATRA-treated vs. untreated mice, at embryonic gestation day 14.5 (E14.5) (n=3 per group). To confirm differentially methylated levels of susceptible genes, real-time quantitative PCR (qPCR) was used to correlate expression of differentially methylated genes related to CP. Results: We identified 196 differentially methylated genes, including 17,298 differentially methylated CCGG sites between ATRA-treated vs. untreated embryonic mouse palatal tissues (P<0.05, log2FC>1). The CP-related genes Fgf16 (P=0.008, log2FC=1.13) and Tbx22 (P=0.011, log2FC=1.64,) were hypermethylated. Analysis of Fgf16 and Tbx22, using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG), identified 3 GO terms and 1 KEGG pathway functionally related to palatal fusion. The qPCR showed that changes in expression level negatively correlated with methylation levels. Conclusions: Taken together, these results suggest that hypermethylation of Fgf16 and Tbx22 is associated with decreased gene expression, which might be responsible for developmental failure of palatal fusion, eventually resulting in the formation of CP.

Abstract: DNA methylation is essential for spatiotemporally-regulated gene expression in embryonic development. TBX22 (Chr X: 107667964-107688978) functioning as a transcriptional repressor affects DNA binding, sumoylation, and transcriptional repression associated with X-linked cleft palate. This study aimed to explore the relationship and potential mechanism between TBX22 exon 5 methylation and palatal shelf fusion induced by all-trans retinoic acid (ATRA). We performed DNA methylation profiling, using MethylRAD-seq, after high throughput sequencing of mouse embryos from control (n=9) and ATRA-treated (to induce cleft palate, n=9) C57BL/6J mice at embryonic gestation days(E) 13.5, 14.5 and 16.5. TBX22 exon 5 was hyper-methylated at the CpG site at E13.5 (P=0.025, log2FC=1.5) and E14.5 (P=0.011, log2FC:1.5) in ATRA-treated, whereas methylation TBX22 exon 5 at the CpG site was not significantly different at E16.5 (P=0.808, log2FC=-0.2) between control and ATRA-treated. MSP results showed a similar trend consistent with the MethylRAD-seq results. qPCR showed the change in TBX22 exon 5 expression level negatively correlated with its TBX22 exon 5 methylation level. These results indicate that changes in TBX22 exon 5 methylation might play an important regulatory role during palatal shelf fusion, and may enlighten the development of novel epigenetic biomarkers in the treatment of CP in the future.